IAU Colloquium 184: AGN Surveys 

ASP Conference Series, Vol. TBD, 2001 

R. F. Green, E. Ye. Khachikian, and D. B. Sanders 



Discovery of Active Galactic Nuclei in Mid- and 
Far-Infrared Deep Surveys with ISO 



Yoshiaki Taniguchi 

Astronomical Institute, Graduate School of Science, Tohoku University, 
Aramaki, Aoba, Sendai 980-8578, Japan 

o. 

CN . Abstract. We present a summary on the discovery of active galactic 

J> \ nuclei in mid- and far-infrared deep surveys with use of the Infrared Space 

O ■ Observatory. 



J> ■ 1. Introduction 

on : 

Active galactic nuclei (AGNs) such as quasars have been one of important issues 
of modern astrophysics since the discovery of a quasar 3C 273 in 1963 (Schmidt 
1963). Their central engine has been thought to be mass-accreting, super- 
massive single black holes around which gravitational energy is transformed into 
huge kinetic and radiation energy with the help of gaseous accretion disks (e.g., 
Rees 1984). Since AGNs are intrinsically bright in any frequency regime, various 
kinds of AGN surveys have been conducted by using radio through optical to 
q X-ray telescope facilities (e.g., Peterson 1997). 

Since current unified models for AGNs introduce a dusty torus around the 
central engine (see for reviews, Antonucci 1993; Urry & Padovani 1995), any 
infrared information is absolutely necessary to understand AGN phenomena be- 
cause such dusty tori radiate anisotropic infrared emission (e.g., Pier & Krolik 
1992; Clavel et al. 2000; Laurent et al. 2000; Murayama, Mouri, & Taniguchi 
2000 and references therein). Another important issue related to infrared ob- 
servations of AGNs seems to be possible starburst-AGN connections (Weedman 
1983; Sanders et al. 1988; Heckman et al. 1989; Mouri & Taniguchi 1992, 2001; 
Cid Fernandes et al. 2001; Storchi-Bergmann et al. 2001). In particular, in- 
frared emission properties are highly useful in investigating circumnuclear star 
formation and its relation to the nuclear activity in AGNs (e.g., Genzel et al. 
1998; Tran et al. 2001). 

Most Seyfert galaxies and a part of quasars were indeed detected at a wave- 
length range between 12 fim and 100 /xm by the Infrared Astronomical Satellite, 
IRAS, launched in 1983 (see for a review, Soifer, Houck, & Neugebauer 1987). 
Although IRAS all-sky survey covered over 96% of the sky, its completeness flux 
limit was 0.5 Jy at 12 fim, 25 /j,m, and 60 /im, and 1.5 Jy at 100 fim. Therefore, 
IRAS could probe infrared galaxies and AGNs up to redshift ~ 0.2 - 0.3 (e.g., 
Veilleux et al. 1999; Borne et al. 2000) except for some unusually bright (or grav- 
itationally amplified) sources beyond redshift z ~ 0.5, e.g., IRAS F10214+4724 
at z = 2.28 (Rowan-Robinson et al. 1991) and IRAS F15307+3252 at z = 0.93 
(Cutri et al. 1994); see also Rowan-Robinson (2000). 
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More than 10 years after the launch of IRAS, the Infrared Space Observa- 
tory, ISO, was launched in 1994 (Kessler et al. 1996). The two infrared array 
detectors, ISOCAM (Cesarsky et al. 1996) and ISOPHOT (Lemke et al. 1996), 
were used to carry out mid-infrared (MIR) and far-infrared (FIR) deep surveys, 
respectively. In this paper, we will give a summary of the MIR and FIR deep 
surveys with ISO and then describe what kinds of AGNs are newly found in 
these surveys (see for a review, Genzel & Cesarsky 2000). 

2. MIR Deep Surveys with ISO 

ISOCAM has two independent channels containing each a 32x32 pixel detector; 
1) the short wavelength channel (2.5 /im to 5.5 /im) and the long wavelength 
one (4 /im to 18 /im) (Cesarsky et al. 1996). Since the short wavelength can 
be accessible from ground-based telescope facilities, the long wavelength camera 
was used to carry out MIR deep surveys; i) at 7 /im (Taniguchi et al. 1997; Sato 
et al. 1999, 2001; Rowan-Robinson et al. 1997; Aussel et al. 1999; Oliver et al. 
2000; Lemonon et al. 1998; Altieri et al. 1999; Fadda et al. 2000), ii) at 12 /mi 
(Clements et al. 1999), and iii) 15 /im (Rowan- Robinson et al. 1997; Aussel et 
al. 1999; Oliver et al. 2000; Lemonon et al. 1998; Altieri et al. 1999; Fadda et 
al. 2000). A summary of these surveys is given in Table 1; note that the 12 /im 
survey by Clements et al. (1999) is not included in this table. 

Among the above MIR imaging surveys, follow-up spectroscopy has been 
done for five surveys given in Table 2. Eleven AGNs (4 type 1 and 7 type 2 
AGNs) were found among 78 sources in the CFRS-ISOCAM survey (Flores et 
al. 1999). However, only a few AGN were found in the remaining four surveys 
(Aussel et al. 1999; Taniguchi et al. 1997; Taniguchi 1999, 2000; Altieri et al. 
1999; Fadda et al. 2000). In total, only 14 AGNs (7 type Is and 7 type 2s) 
were found among 217 sources in the five surveys. This gives an AGN fraction, 
/agn = A^AGN/A^gai ~ 7%. In Figure 1, we show an optical spectrum of a quasar 
at z = 1.025 found by Taniguchi et al. (1997) as an example. 

Although the number of AGNs found in the MIR surveys with ISO is not so 
large, it seems interesting to compare the above observational result with some 
model predictions. Oliver et al. (1997) estimated expected numbers of AGNs 
for their MIR survey of HDF (see also Rowan- Robinson et al. 1997; Aussel et al. 
1999), adopting the following two models; 1) PRR models (Pearson & Rowan- 
Robinson 1996), and 2) AF models (Franceschini et al. 1994); in this article, we 
do not give details of both PRR and AF models. 

a) AGNs expected in the 7 /im survey: In the case of PRR models, the expected 

number of AGNs in the ISO/HDF field is 1.26 (0.66 type 1 and 0.6 type 
2 AGNs) if the 7 /mi limiting flux is 38.6 /iJy and the AGN fraction is 
/agn ~ 22%. On the other hand, in the case of AF models, they are 0.34 
and 9%, respectively. 

b) AGNs expected in the 15 /im survey: In the case of PRR models, the 

expected number of AGNs in the ISO/HDF field is 1.51 (0.73 type 1 and 
0.78 type 2 AGNs) if the 15 /im limiting flux is 255 /iJy and the AGN 
fraction is /agn ~ 21%. On the other hand, in the case of AF models, 
they are 0.25 and 3.6%, respectively. 
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Table 1. A summary of the MIR surveys with ISO 



Field" 


Area 6 


F hm (7u.m) c 


N(7fim) d 


F Um (15/im) e 


7V(15/xm)' 




(arcmin 2 ) 










HDF-N 1 


5/9 


65 


7 


200 


45 


LHNW 2 


9/- 


32 


27 






SSA 13 3 


16/- 


6 


65 






CFRS 4 


-/100 






250 


78 


ELAIS 5 


11.7/22.8 9 


1000 


~ 700 


2000 


~ 800 


A2390C 6 


5.76/5.76 


65 


4 


65 


4 


A2390 7 


6.76/6.76 


25 


31 


40 


34 


A1689 8 


36/36 


150 


41 


300 


18 



a HDF-N = Hubble Deep Field-North, LHNW = the NW field in the Lockman Hole, SSA 13 = 
Hawaii Small Selected Area No. 13, CFRS = Canada-France Redshift Survey Field, ELAIS = 
European Large Area ISO Survey field, A2390C = the central region of Abell 2390, A2390 = 
Abell 2390, & A1689 = Abell 1689. Their references are; 1. Aussel et al. 1999, 2. Taniguchi 
et al. 1997, 3. Sato et al. 1999, 2001, 4. Flores et al. 1999, 5. Oliver et al. 2000, & Serjeant 
et al. 2000, 6. Lemonon et al. 1998, 7. Altieri et al. 1999, & Fadda et al. 2000. 

6 Sky coverage. The first and second numbers means the sky coverage at 7/im and 15 /im, 
respectively. 

c The 3 <r rma flux limit at 7 fim in the survey. 

d The number of objects detected at 7 /im in the survey. 

e The 3 <r rms flux limit at 15 fxm in the survey. 

■^The number of objects detected at 15 fim in the survey. 

9 Note that the area for ELAIS is given in units of square degree. 



Table 2. 


AGNs found in the MIR 


surveys 


with ISO 


Field a 




Aagn c 


/agn 


d A^Typcl 6 A r Type2 / 








(%) 






HDF-N 


49 


1 


2 


1 





LHNW 


13 


1 


8 


1 





CFRS 


78 


11 


14 


4 


7 


A2390 


32 














A1689 


45 


1 


2 


1 





Total 


217 


14 


7 


7 


7 



"The same as those in Table 1. 

6 The number of galaxies detected in the survey. 

c The number of AGNs detected in the survey. 

d The fraction of AGNs = AT AGN /JV ga i. 

e The number of Type 1 AGNs detected in the survey. 

■^The number of Type 2 AGNs detected in the survey. 
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It seems difficult to draw any firm conclusions from the above comparisons be- 
tween the observations and models because of small-number statistics, although 
we do not see significant inconsistency between them. Future MIR surveys will 
give us firmer answers. 

Among the ISO MIR deep surveys, it has been found that there are a 
number of MIR (either 7 ^m or 15 fim) sources without optical/NIR counterparts 
(Taniguchi et al. 1997; Aussel et al. 1999; Flores et al. 1999). These results are 
summarized in Table 3. It is shown that ~ 18% of the MIR sources found in the 
three MIR surveys have no optical/NIR counterparts and ss 10% of them above 
5<7 detection have no counterparts. 

We cannot exclude a possibility that some of such sources may be attributed 
to unexpected noises. However, if a starburst galaxy with a mass of 10 11 M at 
z ~ 3 is heavily reddened (e.g., Ay ~ 10), this galaxy could be detected at MIR 
but not at 2 /im (Taniguchi 2000). Recent submillimeter deep surveys have been 
finding such dust-enshrouded high-z galaxies (Hughes et al 1998; Barger et al. 
1998; Smail et al. 1999). Therefore, we cannot also exclude a possibility that 
deep NIR surveys may miss a certain part of reddened populations, causing the 
underestimate of either star formation density or nonthermal energy density or 
both in the universe. It seems important to keep the presence of such MIR 
sources without counterparts in mind in future investigations. 



Table 3. MIR sources without optical/NIR counterparts 



Field" 


S/N 




N c 

iv no 


f d 

J no 








(%) 


HDF-N 


> 5(7 


49 


7 


14 




> 4(7 


51 


16 


31 


LHNW 


> 5cr 


13 


2 


15 


CFRS 


> 4a 


40 


4 


10 




> 3a 


34 


4 


15 


Total 






187 


34 18 



"The same as those in Table 1. 

6 The number of galaxies detected in the survey. 

c The number of galaxies without optical/NIR counterparts. 

d The fraction of no-counterparts = N no /N ga i. 
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Figure 1 



Figure 1. Optical spectra of four 7 fim sources found in the MIR 
deep survey by Taniguchi et al. (1997). The spectroscopic observations 
were made with use of LRIS on the W. M. Keck 1 telescope (Cowie et 
al. in preparation). One quasar at z = 1.025 was found in their survey 
(see also Taniguchi 1999). The middle panel shows the 7 /im image 
of the LHNW field (halftone; darker is brighter). The contours show 
NIR (HK) image taken with the University of Hawaii 2.2 m telescope. 
Although the seeing size in the NIR image is good (FWHM~0.8 arcsec), 
the image is blurred in order to make the comparison between the 7 
fim and NIR images easier. 
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3. FIR Deep Surveys with ISO 

ISOPHOT is an imaging photopolarimeter covering a wavelength range between 
2.5 /jm and 240 fim (Lemke et al. 1996). Among the several detectors, two de- 
tectors, C100 and C200, were used to carry out FIR deep surveys; i) FIRINDIV- 
DEEP&DEEPf] (Kawara et al. 1998; Matsuhara et al. 2000), ii) FIRBACK- 
ELAIS; the 90 /xm survey (Efstathiou et al. 2000; Serjeant et al. 2001), and the 
170 fim survey (Puget et al. 1999; Lagache & Dole 2001; Dole et al. 2001). A 
summary of these survey is given in Table 4; note that we do not include another 
ISOPHOT survey at 60/um and 90 p by Linden- Vornle et al. (2000) because 
no follow-up observation is available. 

Although optical followup spectroscopy has not yet been completed for the 
two surveys, some preliminary results were reported recently. Serjeant et al. 
(2001) made optical spectroscopy for 20 sources detected in the FIRBACK- 
ELAIS survey with /(90/um) > 100 mJy, which are also detected either at 15 
/j,m or at 1.4 GHz. Among them, they found two Seyfert galaxies at z = 0.149 
and z = 0.225; the remaining 18 sources are either starbursts (16 sources) or 
early- type galaxies (2 sources). This gives an AGN fraction, /agn = 2/20 = 
10%. On one hand, Kakazu et al. (2001a, 2001b; see also Murayama et al. 2001) 
made optical identification of 35 sources found in the FIRINDIV-DEEP&DEEP 
survey with /(170/xm) > 100 mJy, using the W. M. Keck, Subaru, UH88, and 
VLA facilities. They identified two AGNs; a Seyfert 1.5 galaxy at z = 0.206 
and a quasar at z = 1.60. The remaining 33 sources are classified as starbursts 
(21 sources), LINERs (10 sources), and early-type galaxies (2 sources). Since 
infrared-selected LINERs may be shock-heated galaxies (Taniguchi et al. 1999; 
Lutz, Veilleux, & Genzel 1999; cf. Imanishi, Dudley, & Maloney 2001), they are 
not genuine AGNs but starburst-related (i.e., superwind) galaxies. Therefore, 
their survey gives an AGN fraction, /agn = 2/36 6%. 

The above AGN fractions are considered to be tentative values because the 
number of observed galaxies are still small. We hope that future followup obser- 
vations will give us firmer statistics. It seems worthwhile noting that Frances- 
chini et al. (1989) estimated an AGN fraction in such FIR deep surveys; i.e., 
/agn ~ 10%, being similar to the observed values. 

Finally, we mention that 10 ultraluminous infrared galaxies (ULIGs) and 1 
hyperluminous infrared galaxy (HyLIG) were identified as counterparts of 170 
//m sources found in Kawara et al.'s (1998) survey (Kakazu et al. 2001a, 2001b; 
Murayama et al. 2001), Although the HyLIGg is a quasar at z = 1.6, the ULIGs 
are located at z ~ 0.3 - 0.8. Such intermediate- z ULIGs have not yet been found 
by IRAS. It is also remarkable that the observed surface density of ULIGs is ~ 
10 degree -2 . Note that only ~ 120 ULIGs with z < 0.3 were found by IRAS 



1 Contrasting to the FIRBACK-ELAIS program (acronym for Far Infrared background), 
in this article, we call the Japanese-IfA/UH program FIRINDIV-DEEP&DEEP where 
INDIV means individual sources found in the FIR survey and DEEP&DEEP means a 
deep survey for dust- enshrouded extragalactic populations (deep). Please make sure that 
DEEP&DEEP is different from the famous DEEP (Deep Extragalactic Evolutionary Probe) 
survey program promoted by David Koo at UCO/Lick Observatory, University of California 
( http: / /www / ucolick.org~deep /home. html ) . 

2 See for discovery of a HyLIG in the ELAIS survey, Morel et al. (2001). 
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(Veilleux et al. 1999; Borne et al. 2000; see for a review Sanders &: Mirabel 
1996), giving a surface density of ~ 0.003 degree -2 . Therefore, it is suggested 
that the surface density of intermediate- z ULIGs is much higher by a factor of 
3000 than that of low-z ULIGs. Since it has been often argued that ULIGs 
may be precursors of quasars (Sanders et al. 1988; Sanders & Mirabel 1996; 
Taniguchi, Ikeuchi, &: Shioya 1999), it is interesting to investigate the nature 
of intermediate- z ULIGs found in the ISO FIR deep surveys. Future wide- field 
FIR deep surveys and their follow-up observations will be also very important 
to find new populations of intermediate- z and high-z ULIGs and HyLIGs. 



Table 4. A summary of the FIR surveys with ISO 



Program 


Field" 


F lim (90^m) b 
(mJy) 


iV(90/im) c 


F lim (170 / um) d 
(mJy) 


N(170u,m) e 


FIRINDIV- 


LH 


45 


36 


45 


45 


DEEP&DEEP 


1.1 deg 2 / 
1.1 deg 2 










FIRBACK- 


Several 


100(?) 


120(?) 


135 


106 


ELAIS 


11.6 deg 2 / 
3.89 deg 2 











a LH _ Lockman Hi Hole. As for the FIRBACK survey fields, see http://wwwfirback.ias. u 



psud.fi. The first sky coverage is for the 90 fim survey and the second one is for the 170 fxxa 



one. 
6 3a 

rms - 

c The number of galaxies detected in the survey. The galaxies detected in FIRINDIV have 
F(90/mi)> 150 mJy. 

3<T rms • 

e The number of galaxies detected in the survey. The numbers given in this column are that of 
galaxies with _F(170/^m)> 150 mJy for FIRINDIV and that of galaxies with F(170/im)> 180 
mJy for FIRBACK. 



4. Discussion and Summary 

As summarized in this article, ISO enabled us to perform a number of MIR and 
FIR deep surveys (see also Genzel &: Cesarsky 2000). We give a summary of 
main points of this article below. 

1) As for the MIR surveys, the ISO deep surveys at 7 //m and 15 fim are 
sensitive enough to detect sources down to 10 fiJy (Taniguchi et al. 1997; 
Aussel et al. 1999; Sato et al. 2001b). Approximately, 10% of the detected 
sources appear AGNs from low-z Seyfert galaxies to a high-z quasar at z 
1-2 (Taniguchi 1999; Aussel 1999; Flores et al. 1999). Since sky areas 
observed by the surveys are so small, it seems hard to make statistical 
arguments on AGN populations. However, number counts of the MIR 
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sources show significantly stronger evolution than what expected from no 
evolution models (e.g., Elbaz et al. 1999). This may be attributed mainly 
to intense star formation in galaxies at intermediate- and/or high-z galax- 
ies (e.g., Takeuchi et al. 2000; Chary & Elbaz 2001). Since dusty starburst 
galaxies sometimes harbor hidden AGNs (e.g., Sanders et al. 1988; Ivison 
et al. 2000; Willott et al. 2001), future follow-up observations will be 
important to understand what faint MIR sources are. 

2) The ISO FIR deep surveys have also shown that approximately, 10% of the 

detected sources appear AGNs (Serjeant et al. 2001; Kakazu et al. 2001a, 
2001b; Murayama et al. 2001). Although the total sky area surveyed 
by the FIR surveys exceeds 10 deg 2 , optical follow-up observations have 
not yet been fully done. Therefore, future follow-up observations will be 
important to understand the nature of faint FIR sources. The most re- 
markable finding of the ISO FIR deep surveys is the discovery of numerous 
ULIGs at intermediate redshift between z ~ 0.3 and z ~ 0.8 (Kakazu et 
al. 2001a) because such populations have not yet known from the IRAS 
survey. It seems quite likely that these populations contribute to the ob- 
served excess number count at FIR (e.g., Elbaz et al. 1999; Takeuchi et 
al. 2000). 

3) The above exciting results urge us to conduct new MIR and FIR surveys. In- 

deed, new space infrared telescope facilities will be launched soon; SIRTF| 
and IRIS (ASTRO-FQ). We hope that these facilities will give us much more 
information on AGN populations in the infrared universe. 
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